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Perturbed  Bifurcation  of  Stationary  Striations 
in  a  Contaminated,  Non-uniform  Plasma 


J.  F.  Magnan  and  R.  A.  Goldstein 

Abstract 

A  cylindrical,  weakly  ionized,  and  collision  dominated  neon  plasma 
described  by  a  system  of  nonlinear,  parabolic  reaction-diffusion  equations 
for  the  electron  and  metastable  atom  axial  densities  exhibits  a  bifurcation 
from  a  uniform  to  a  striated  state  at  a  critical  length  of  the  plasma 
column.  The  sharp  transition  between  states  predicted  by  the  theory  is 
in  contrast  with  the  smooth  transition  observed  in  experiments.  We  apply 
the  theory  of  singular  perturbations  of  bifurcations  to  show  that  small 
inhomogeneities  in  the  plasma,  such  as  those  caused  by  non-uniform  heating 
and  contamination,  are  sufficient  to  qualitatively  explain  the  experimental 
results.  Ue  observe  that  a  steady,  axial  magnetic  field  in  the  plasma 
can  also  produce  a  smooth  transition. 
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1.  Introduction 

It  has  been  found  that  the  longitudinally  uniform  and  axially  symmetric 
steady  state  of  a  cylindrical  plasma  may  bifurcate  transcritically  to  a 
non-uniform  striated  state  after  the  onset  of  an  ionization  instability  [l]. 
The  weakly  ionized  and  collision  dominated  noble  gas  plasma  which  is  uniformly 
heated  with  radio- frequency  electromagnetic  waves  can  be  described  by  a 
system  of  non-linear  reaction-diffusion  equations  for  the  electron  N  and 
metastable  atom  M  population  densities.  The  theory  in  [l]  consider  the 
following  initial -boundary  value  problem  for  N  and  M: 


(1.1a) 

Nt  “  DNzz  +  F(N,M) 

,  0  £  z  £ 

L 

(1.1b) 

M  -  9 DM  +  G(N,M) 

c  zz 

,  0  £  z  * 

L 

(1.1c) 

N  -  No  .  M  -  Mo 

at  z  *  0  and 

z  *  L  , 

(l.ld) 

N(z,0)  -  h(z,e )  , 

M(z,0)  -  i(z,e) 

• 

In  (1.1)  t  is  the  time  variable  and  z  the  axial  coordinate  in  a  cylinder 
of  length  L.  D  is  the  ambipolar  diffusion  coefficient  of  the  ions  and  0 
is  the  ratio  of  the  metastable  atom  diffusion  coefficient  to  D.  and  Mq 
are  the  uniform  steady  state  values  of  N  and  M.  The  initial  functions  h 
and  i  have  asymptotic  representations  eh^(z)  and  ei^(z)  which  have  Fourier 
sine  series  expansions  and  vanish  at  the  boundaries.  The  small  parameter  s 
is  defined  below.  The  functions  F  and  G  are  nonlinear  particle  production 
and  loss  rates  describing  the  interactions  between  electrons,  metastables, 
and  neutral  gas  attxas  in  the  plasma.  They  are  discussed  in  [l],  [2],  [3], 


Their  form  is  given  in  [2]  as 
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(l.le)  F(N,M)  =■  -  f  K  +  f2NM  +  , 

(l.lf)  G(N,M)  -  gx  -  g2M  -  g3NM 

where  the  positive  constants  f^,  f2,  f^,  g^,  g2,  and  g^  are  rate  coefficients 
obtained  by  properly  averaging  collision  cross  sections  over  the  electron 
energy  distribution  function  in  the  plasma. 

An  asymptotic  expansion  of  the  solution  to  problem  (1.1)  as  e  -*  0  is 
obtained  in  [l]  by  a  standard  two-time  method.  The  analysis  of  the  problem 

finds  that  the  uniform  state  may  become  unstable  when  a  critical  value  u 

L2  C 

of  the  bifurcation  parameter  y  =  — r  is  exceeded.  For  y  near  y  an  initial 

Dtt  c 

perturbation  of  the  uniform  state  will  develop  into  a  non-uniform  steady 
state  with  a  sinusoidal  density  variation  along  the  plasma  axis.  The 
striated  state  of  the  plasma  is  given  by 


(1.2a) 


"n  " 

a 

i 

as 

o 

_ 1 

+  eA 

*i  - 

_  M  _ 

o* 

c 

A  . 

c 

ncnZ  2 

sin  — —  +  0(e) 


there  the  critical  mode  n 

c 


«n 

(1.2b)  Ac - ~M'(0) 


is  odd, 


9 


9 


and  ,  a,  b^  are  constants  defined  below.  The  normalization  y  -  y^  -  c 

c  nc  c 

defines  e  and  u'(0).  The  value  of  n  in  (1.2a)  is  determined  by  u  since 

c 

disconnected  regions  of  this  parameter  may  exist  where  modes  of  different 
n  grow  in  amplitude.  For  even  n,  it  is  necessary  to  choose  y' (0)  =  0  and 
a  higher  order  calculation  is  required. 

Experimental  investigations  of  the  phenomenon  have  qualitatively  con¬ 
firmed  the  results  of  the  analysis  [4].  Experiments  [5]  show  however 

that  the  transition  from  the  uniform  to  the  striated  state  (n  -  1)  is 

c 
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smooth,  in  contrast  with  the  theoretical  prediction  that  an  abrupt  change 
of  state  occurs  at  the  bifurcation  point.  The  contrast  between  theory  and 
experiment  is  shown  by  Fig.  1.  We  attribute  this  discrepancy  to  the  non- 
uniform  heating  of  the  plasma  by  the  radio- frequency  electric  field  and  to 
the  presence  of  impurities  in  the  plasma.  In  this  paper  we  show  that  these 
effects  are  sufficient  to  qualitatively  explain  the  smooth  transition 

between  states.  In  addition,  we  show  that  the  transition  can  be  smoothed 
by  a  steady  magnetic  field  in  the  plasma. 

In  §2  we  modify  problem  (1.1)  to  include  the  effect  of  impurities  and 
non-uniform  heating  which  are  assumed  to  be  of  the  same  small  amplitude  6 • 

We  specifically  consider  the  case  of  an  argon  contaminant  in  a  neon  plasma 
and  call  the  resulting  steady  state  problem  the  perturbed  problem. 

We  apply  the  method  of  Matkowsky  and  Reiss  [6],  [7]  in  §§3-5  to  find 
asymptotic  solutions  of  the  perturbed  problem.  The  technique  used  is  that 
of  matched  asymptotic  expansions  where  an  outer  and  inner  boundary  layer  ' 
expansion,  which  are  perturbation  expansions  in  5  of  the  bifurcation 
branches  away  from  and  near  to  the  bifurcation  point,  are  connected  by  a 
matching  procedure.  A  uniform  asymptotic  representation  of  the  solution  can 
then  be  obtained.  The  method  has  been  applied  to  the  buckling  of  columns, 
forced  nonlinear  oscillations,  Rayleigh-Benard  convection,  Poiseuille  flow, 
and  panel  flutter  [7],  [8],  [9],  [10].  In  our  analysis  we  have  retained 
much  of  the  notation  used  in  [l]. 
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2.  The  Perturbed  Problem 

Impurities  in  a  laboratory  plasma  are  always  present.  They  consist 
mainly  of  water  vapor  and  air  from  leaks  and  surface  desorption.  Two  sample 
reactions  with  impurities  are 

(2.1)  e  +  H20  -»  H"  +  OH  , 
and 

(2.2)  A  +  Ne  -»  A+  +  Ne  +  e  , 

in 

where  the  two  processes  are,  respectively,  dissociative  attachment  and 

Penning  ionization  (of  an  argon  ground  state  atom  A  by  a  neon  metastable 

atom  Ne  ).  Each  impurity  has  a  different  effect  on  the  plasma:  water 
m 

molecules  destroy  electrons;  and  argon  atoms  create  electrons  and  destroy 
metastable  atoms. 

We  shall  consider  a  neon  plasma  contaminated  by  an  argon  impurity  which 
is  introduced  at  a  constant  rate  by  an  external  source.  The  argon  reacts 
almost  immediately  with  the  metastables  because  of  the  large  argon  collision 
frequency  and  high  probability  for  Penning  ionization  to  occur.  For  this 
reason  we  can  neglect  the  direct  ionization  of  argon  by  electrons,  the 
diffusion  of  argon,  and  the  change  of  the  argon  reaction  rate  with  small 
variations  in  the  metastable  concentration.  Under  these  conditions  the 
conservation  equation  for  the  argon  concentration  is  decoupled  from  the  N 
and  M  equations  and  gives  the  result  that  the  steady  state  concentration  has 
the  same  spatial  dependence  of  the  source. 

The  time  rate  of  change  of  the  electron  and  metastable  atom  populations 
due  to  an  argon  impurity  concentration  6I(z)  is  proportional  to  SI(z)M, 
where  5  >  0  is  a  small  parameter  characterizing  the  magnitude  of  the  impurities. 


i 
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It  may  be  defined,  for  example,  as 
(2.3)  6  *  max 

In  the  theory  presented  in  [l]  it  is  assumed  that  the  plasma  is  uni¬ 
formly  heated  by  the  radio- frequency  electric  field  so  that  the  reaction 
rates  do  not  vary  along  z.  In  experiments,  however,  the  heating  can  be 
non-uniform  and  the  rate  coefficients  of  the  plasma  depend  on  z.  In  the 
experiments  described  in  [4]  the  heating  is  concentrated  near  the  ends  of 
the  plasma  where  the  field  produced  by  two  end  electrodes  is  strongest. 
The  large  field  at  the  ends  is  due  to  the  proximity  of  the  electrodes  and 
the  shielding  of  the  central  region  of  the  plasma  by  the  plasma  at  the 
boundary. 

If  we  assume  the  non-uniformity  of  the  field  to  be  small  and  of  the 
same  characteristic  magnitude  6  as  the  contamination  in  the  plasma  we  may 
write  the  space  dependent  part  of  the  reaction  rates  as 

(2.4a)  6  [  -cr1(z)N  +  «2(z)NM  +  a3(z)N2M]  =  6X(z,N,M)  , 

and 

(2.4b)  8[01(z)  -  02(z)M  -  03(z)NM]  ■  6Y(z,N,M) 

where  the  a's  and  (3's  are  rate  coefficients. 

The  steady  state  problem  for  the  non-uni formly  heated,  contaminated 
plasma  is  given  by 

(2.5a)  N?s  +  pi  F(N,M)  +  6(kMI(5)  +  X(|,N,M))]  -  0  ,  0  £  l  «s  tt  , 
(2.5b)  GMgg  +  pi  G(N,M)  +  6(-kMI(?)  +  Y(§,N,M))]  -  0  ,  Os^Srt  , 

(2.5c)  N  -  NQ  ,  M  -  MQ  at  ?  -  0  and  ?  -  tt  , 


TTZ 

where  5  *  —  is  a  dimensionless  space  variable  and  k  is  a  rate  coefficient. 
Equations  (2.5)  define  the  perturbed  bifurcation  problem  when  6  j4  0;  and 
the  bifurcation  problem  when  6*0. 
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3.  The  Outer  Expansions 

The  four  sub-branches  of  the  bifurcation  branches  are  shown  in  Pig.  1. 

The  branches  include  the  uniform  state  sub-branches  to  the  right  and  left 

of  and  the  transcritical  striated  state  sub-branches  branching  above  and 
c 

below  u  .  Let  the  vector  u  with  components  N  and  M  be  the  outer  expansion; 

*  c  ^ 

and  let  un  be  a  bifurcation  sub-branch.  We  seek  asymptotic  expansions  as 
6  -*  0  of  the  solution  of  the  perturbed  problem  near  in  the  form 


GO 

(3.1)  u  -  uj(§,M)5J 


The  coefficients  u.  of  the  outer  expansion  are  obtained  by  substituting 
(3.1)  into  (2.5)  and  equating  to  zero  coefficients  of  6^.  This  results  in 
the  following  linear  boundary -value  problem  for  u^: 

(3.2a)  (L  +  /lOu^  ■  R  ,  u^  *  0  for  §  =  0  and  §  »  rr  , 

where 


(3.2b)  L 


F  F 
N  M 


G  G 
N  M 


(3.2c)  R  -  (5 ) 


”-l" 
.  1  . 


X(§  ,Nq,Mq) 

y(5,n«,m'). 


In  (3.2)  the  derivatives  of  F  and  G  are  evaluated  at  and  N^,  are 
the  components  of  u^. 

We  solve  (3.2)  by  using  eigenfunction  expansions.  The  eigenfunctions 
vQ  and  eigenvalues  of  the  operator  L  +  fiQ  are  defined  by 
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(3.3)  (L  +  jiQ)v  *  a  v  ,  v  -  0  for  |  =  0  and  . 

^  n  n  n  n  1  r 

To  solve  the  eigenvalue  problem  (3.3)  for  the  uniform  branch  where 
Nq  *  Nq  and  we  rewrite  (3.3)  as 

(3.4)  [L  +  McQq  +  <H  -Mc)Q0>n  *  <Vn  , 

where  QQ  is  Q  evaluated  at  N^,  M^.  The  eigenfunctions  v  ,  v  and  eigenvalues 

^  —  ll  ^ 

e^,  e^  of  the  operator  L  +  pQQ  and  its  adjoint  L  +  jiqq  are  defined  by 

(3.5a)  (L  +  /jQq)v±  *  e^vn  ,  v”  =  0  for  §  *  0  and  §  =  tt  , 
and 

(3.5b)  (L  +  =  e^Vj^  ,  v^  =»  0  for  §  -  0  and  5  »  tt  , 


where  e  is  the  complex  conjugate  of  e  .  In  (3.5)  each  operator  has  two 
n  n 

families  (+  and  -)  of  eigenfunctions  and  eigenvalues.  They  are  given  in 
Cl]  as 


(3.6a)  v+ 
n 

(3.6b)  v“ 
n 


sin  ng 


sin  n§ 


bn  ‘  (en  +  "2  '  “V^M  • 
cn  '  (en  +  1,2  '  ■ 


(3.6c)  v*+  ■  E 


sin  n?  ,  E  =  - - - - * 

n  "l  +  bb 
n  n 


bn  '  +  1,2  '  > 


(3.6d)  v  ”  “  K 
n  n 


1  “i 
* 

Lcn  J 


sin  ng  ,  K  »  J- - - - r 

n  n  1  +  c  c 
n  n 


cn  *  <*o  +  "  -  , 


9 


9 


HP 
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(3.6e)  e;.|  ,![!„*  (t2-4P„)!e]  , 


n  n  n 
2 


<3.6£)  In  -  ^  +  F„  -  J-  (1  +  6)  ,  Pn  -  *J-  -  (8F„  +  V  f  +  -  FmGn 


In  Cl]  it  was  assumed  that  T  <  0  and  P  ( — •)  =  0. 

n  n\u  / 

'c 

The  eigenfunctions  satisfy  the  orthogonality  conditions 


(3.7a)  <v\  v*>  =  0  for  n  +  m  , 

n  m 


(3.7b)  <v*  v*)  -  0  , 

n  m 


(3.7c)  (v^.v*)  -  1  , 

'  n  n  5 


where  the  inner  product  of  two  vectors  a ,  b  is  defined  by 

.  <"  .t, 


(3.7d)  (a,b>  »  (  J(5)b(5)<i5 


and  a^  is  the  Hermetian  adjoint  of  a. 

We  denote  the  eigenfunctions  of  the  operator  L  +(i  Q.  and  its  adjoint 


cx0 

'k  i  ^ 

L  +  McQq  by  Xn  an<*  respectively. 

i  I 

It  is  useful  for  the  matching  in  §5  to  expand  the  v  and  v  *  in  terms 

n  n 

of  the  v  and  v  as 
~n  ~n 


(3.8a)  V*  -  f  «&£  +  8^)  . 

k=l 


(3.8b)  v  -  r  (a  ^  +  8^  )  , 

k*l 


where 


■  <:£■%>  • 


(3.8c)  <4  -  <£,v*> 


> 
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and 


(3.8d)  a 


*± 

'nk 


*+  *±. 
<v  ,v  > 
Me  n 


3*? 

nk 


<2k  -v„  > 


The  expansion  coefficients  in  (3.8)  are  functions  of  p  and  are  given  by 


(3.9a)  a+  =  a+ 


1  +  bn(M)bn0ic) 

“  1  +  bn(Mc)b>c) 


a  s  a 
n  nn 


1  *  Cn^)bn^c> 

1  +  WW 


(3.9b)  p+  s  p+ 


1  +  bn^i)cn^c) 

0X1  1  +  c  Qi  )c*(u  ) 

n^c  n^c 


8  ■  p 

n  nn 


1  +  ca»>cn<*te) 
1  +  Cn(^c)cn(^c) 


(3.9c) 


*4.  *4-  1  +  b  (^  )b  fyt) 

*+  _  *+  n^c  n^ 

a_  =  a =  — - —  1  - 

1  +  bn^)bn^> 


nn 


*-  1  +  K&J'nW 

a  =  cc  «*  . .  ■—  . 

n  nn  .  .  ,  .  *,  . 

1  +  cnQi)cnQi) 


*+  *+  1  + 

(3.9d)  B  +  M,! - S . C-T* - 

1  +  bn(M>bn(M> 


*_  *_ 

8  3  8 

n  nn 


1  +  CnOOc..(H) 

_ n  c  n 

1  +  cn(^)c*(^t) 


and 


(3.9e)  a 


nk 


*±  ±  *± 

a, =8,  =  8  ,  =0  for  k  ^  n 
nk  nk  nk 


Thus  from  (3.8)  and  (3.9) 


(3.10a)  v*  -  a*v\  +  8*v 
n  rwi  n~n 


and 


(3.10b)  v 


*±  *+  *±  *- 
■  a  v  +  8  v 
n  ~n  n  ~n 


The  asymptotic  behavior  of  the  coefficients  as  p  -*  p  is  needed  below 
and  is  given  by 


(3.11a)  cT  -  8_ 
n  n 


C  -  P*"  ■  1  +  0Qx  -  a  )  » 


(3.11b)  a.n  -  Pn  =  an  =  Pn  =OQi-uc) 

Now  (3.4)  may  be  written  as 

(3.u)  cV  -  +  ^  ^ +  ■ 

where  we  have  used  (3.10a)  and  (3.5a).  Taking  the  inner  product  of  (3.12) 
with  and  using  the  orthogonality  of  the  eigenfunctions  we  arrive  at  an 
expression  for  the 

(3.13a)  c£  -  » 

where 

(3.13b)  £  »  vre>c)  +  9nCe>C)  * 

(3.13c)  <£  -  aV*q*  +  aW  +  a*Vy+  +  3*0  *V 

w  ;  'n  n  n  ^  nyn  4n  n  n4n  n  4n  ’ 


(3.13d)  <  •  <vv£>  > 


.  ±  Tv 

<VVn> 


We  expand  the  inhomogeneous  term  R  in  the  v~  as 


•U)  ~  ‘  j?i(rjvj  +  Vi 


)  . 


The  expansion  coefficients  r  are  given  by 


(3.15a)  r*  -  <v*,R> 


and  have  the  values 


(3.15b)  r*  -  atV  +  P^rT  , 


where 

(3.15c)  Rj 
(3.15d)  R” 

(3.15e)  B 

(3.15f)  Ij 
(3.15g)  Sj 

(3. 15h)  Tj 


ff  "  Sj)  • 

|»‘CJ(».0IJ-TJ)  , 

‘’tO'c1-1  _  c 

l  +  ’  i  1  +  Cj(Mc)cj(Mc) 

^I^sin  jgdg  , 


-  Uxn<5)  +  b?0O  V^sin  JS  dg  . 

bj(Mc)-i3o  j 


iH1 -  <V«)  +  c*Qi  )Y0(S)]8in  jg  dg  . 

Cj<Mc)-l]0  J  C 


and 


9 


(3.15i)  XQ(g)  s  X(g,N0,M0)  ,  YQ(g)  3  Y(g,N0,M0) 


We  assume  R.  £  0. 
■^c 


The  solution  of  (3.2)  is  given  by 


(3.16) 


±  ±  ± 

where  the  Vj,  <Jy  and  Tj  are  given  in  (3.10a),  (3.13a),  and  (3.15b). 
outer  expansions  corresponding  to  the  uniform  sub-branches  are 


(3.17)  u 


N„ 


<-M0 


+  5u.  +  0(6  )  , 

^1 


The 


where  u.  is  given  by  (3.16). 
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The  asymptotic  behavior  of  (3.17)  as  p  -*  p  is  given  by 
(3.18)  u  - 


"V 

+  1u&  f  > 

Bi(kM0Ii-Si>  +0<M-Mr) 

r*  l 

-Mo  - 

e;+(,-Mc)q;+e;°(M-Mc)+oc  ^-mc>2] 

/■N 

o 

3 

X3 

_ J 

ci(kMoW  + 


ej  +  +  ej0(^  ”MC)  +  °C  (M-Mc>  J 

0[  6Qi-«  )]  +  0(52)  . 


VMc>J 


sin  j§  + 


For  j  ■  j  the  first  term  inside  the  brackets  in  (3.18)  becomes 


(3.19) 


BJ  <kVi  -sj  >  r 

c  Jc  ~*c 


sin  jc§  +  0(1) 


since  e  ■  0.  It  is  seen  from  (3.19)  that  the  outer  expansions  (3.17) 
Jc 

are  unbounded  at  the  bifurcation  point  \i  «*  fi  .  They  are  shown  in  Fig.  2 

c 

for  aR.  >0  and  aR+  <  0. 

Jc  Jc 


The  outer  expansions  for  the  striation  sub-branches  (1.2)  will  depend 
on  e  because  of  the  dependence  u^,  /i  and  Q  now  have  on  e .  The  eigenvalue 
problem  (3.3)  thus  becomes 


(3.20)  [L  +  ^(e)Q(e) ]v  (e)  =»  a  (e)v  (e)  ,  v  =0  for  5=0  and  §  -  n  . 

n  n  n  n 

The  matrix  Q  is  evaluated  at  the  striation  sub-branch  and  has  the  form 
(3.21a)  Q(e)  ■  QQ  +  cA  D  sin n  {  +  0(e2)  , 


where 


4 


(3.21b)  D 


lGm  +  b  C 


JNN  +  bn  GNM  GNM 
c 


We  solve  (3.20)  by  seeking  perturbation  expansions  in  e  of  v  and  a  of 

n  n 

the  form 


CD 

1.22a)  v  (e)  *  J'  v  ,c 
j-0  n'j 


(3.22b)  a  (e)  -  £  cr  .ej 
n  j=o  n‘j 

The  coefficients  v  ,  a  in  (3.22)  are  obtained  by  substituting 

*  J  n*  J 

(3.22)  into  (3.20)  and  equating  to  zero  the  coefficients  of  e^.  This  leads 

to  the  following  equations  for  v  .  and  v  - : 

n,0  ~n, 1 

(3.23)  (L  +  a  Q_)v  -  ■  a  _v  _ 

*V*0  n,0  n,0  n,0 


(3.24)  a+McQ0-»ni0)Vn>1  -  -McAcD.i»nc5)viii0 

Equation  (3.23)  is  Che  same  as  (3.5a)  and  therefore  v  „  »  v± .  and 

n,0  ~n’ 

s  n  m  e„*  Thus  we  can  rewrite  (3.24)  as 
n,u  n 

(3.25a)  (L  +  McQq  "  en^Vn  1  *  ^ sin  n?  +  z  si-nnc?  sin  n?  , 
where 


(3.25b) 


ffn,l  -M'(0)(FN  +  bnFm) 


Vn.l  '  +  b„V 


15- 


and 


(3.25c) 


“  Z- 

"  F _ +  b  +  b  F^, " 

zi 

NN  n  NM  n  NM 

Z  3 

■  ■  (J  A 

c 

^c  c 

G _ +  b  G.„,  +  b  G.,., 

Lz2  ^ 

L  NN  n  NM  n  NM  J 

c 


It  will  not  be  necessary  to  find  the  v  ,  in  (3.25a)  because  they 

n,i 

ultimately  give  terms  of  0(e£)  in  the  outer  expansions  and  are  of  higher 

order  than  the  0(c)  and  0(5)  terms  we  are  keeping.  The  a*  1  are  obtained 

n,  i 

from  the  solvability  requirement 


(3.26)  (Attain  ng  +  Zsinn  S  sin  n§)  *  0 

Ml  C 

This  condition  gives  us  equations  which  may  be  solved  algebraically  for 

.  ± 
the  a 

n,l 


(3.27.)  f  („;  +  b*n„;>  + 1  (J-  +  -  &r+v)  <n +  b>2>  - 0  • 

c  c  c 

(3.27b)  &  (w7  +  cV)  +  f  (—  +  - -  -  T-~ - )  (z.  +  c*z_)  =  0  , 

21  n  2  4  \n  2n-n  2n  +  n  /  1  n  2 


for  even  n  and  odd  n  j*  n  , 

c 

and 


(3.27c>  B  (»t  *  b*«£  +  ^  (*!  +  »*z2)  -  0  , 


(3.27d)  B  O'  +  c*w2)  +  V2>  -  0  , 


for  n  ■  n  .  We  solve  (3.21c)  for  a  .  (it  will  be  needed  later)  and  obtain 
c  nc,l 

(3.28)  .--m'<0)<£ 

c’1  c 
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The  solution  of  (3.2)  for  the  striatlon  sub-branches  is  given  by 

'at  .  R' 


(3. 


-1  3 


ijj  i  i 

where  the  v.,  <3  ,  and  R.  are  given  by  (3.22)  and  (3.15c,d).  Thus  it  follows 
J  J  3 

from  (3.1)  and  (3.29)  that  the  outer  expansions  corresponding  to  the 
striation  sub-branches  are 


(3.30)  u  = 


r  n„  i 

"  1 

0 

1  +  eA 

1  M  c 

b 

L  0  J 

n  J 

U  c  _J 

sin  n  §  + 
c 


i*  I 

j=l 


TT 


w.-v 

’  1  " 

C.ttM,!,-!,) 

'I  ' 

1 

+  ^  + 

-bj- 

•I  *  -1.1 

-Cj- 

4 

sin  j§  + 


0(«6 )  +  0(e2)  +  0(62) 


At  the  bifurcation  point,  where  e.  “0  and  e  *  0,  the  outer  expansions 

Jc  +  + 

(3.30)  are  unbounded.  They  are  shown  in  Fig.  2  for  crR^  >  0  and  arR  <  0. 


We  have  used  (3.28)  to  draw  the  curves 
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4.  The  Inner  Expansions 

These  expansions  are  valid  near  p  as  6  -»  0  and  are  obcained  by  stretching 
the  neighborhoods  of  jz  with  the  transformation 

®  . 

(4.1)  nM  -  u  +  H**  +  Z  VV“)J 
c  j=2  J 

where  1]  is  the  stretched  variable  of  the  method  of  matched  asymptotic 
expansions  and  the  small  parameter  v  is  defined  by 


(4.2)  6 (v)  -  vb  . 

The  constants  a  and  b  are  positive  integers  since  we  require  that  fz(v)  -* 
and  6  (v)  ■+  0  as  v  -*  0  and  that  the  derivatives  of  /z(v)  and  6(v)  be  bounded 
as  v  -*  0. 

Let  the  vector  iu  with  components  n(v)  «  N(5  ,/z(v)  ,5  (v))  and 
m(v)  «  M(5,jl(v),6  (v))  be  the  inner  expansion.  We  seek  asymptotic  expansions 
of  the  solution  of  the  perturbed  problem  near  the  singular  point  defined  by 

(4.3)  M  31  Mc  a**  “fo  * 
in  the  form 

0®  . 

(4.4)  U)  -  y  III ,(S,M(v),S(v))vJ 

~  j-0 


where  ui.  has  components  n.,m,. 

~J  J  J 

In  the  usual  manner  we  arrive  at  the  following  equations  for  the  w. : 

/  %  J 

a  £i 

LV5'-!  ] 

ui.  »  0  for  5*0  and  5  ■  tt  , 


r  n 

f-1  1 

x0(5) 

*Mc6'  ■ 

kM  K5) 

u  l  1  J 

/■s* 

O 

>« 

* 

L 


(4.6)  (L  +  HcQ  0)(U2 


+ 


-  M’QoiSi  -  Mc 


2  FNN  +  nlmlFNM 


2  GNN  +  nlmlGNM 


Mc6' 


kr^Kg) 


--ii  rx^sn 

.  ij  L V5>J 
1 


r-ii 

”X0®  >" 

■)  ■ 

kM_I(g) 

°  L  lj 

.Y0(g). 

*  fe 


(u-  -  0  for  §  *  0  and  §  ■  tt  , 


and  in  general  for  j  >  2, 

(4.7)  (L  +  /icQ0)(Uj  ■  ,  a) j  -  0  for  g  »  0  and  g  *  tt 

The  primes  on  6  and  ji  denote  derivatives  with  respect  to  v  evaluated  at 
v  •  0  and  come  from  Taylor  expansions  of  6  and  fJ,.  In  (4.6)  X^(g)  s  X(g,n^,m^) 
and  Y1(g)  s  Y(g,n1,m1). 

The  operator  L  +  not  invertible  and  thus  the  inhomogeneous 

problems  (4.5)-(4.7)  have  solutions  if  and  only  if  jj^  (j  =  1,2,* ••)  is 

•fg  *fg  | 

orthogonal  to  the  null  space  of  L  +  fJcQg  which  is  spanned  by  v^  .  The 

c 

solvability  conditions  for  these  problems  are  thus 


(4.8)  <5^  ,£.)  -  0  ,  J  -  1,2,— 

c  J 

For  J  -  1,  (4.5)  and  (4.8)  imply  that 

(4.9)  6'B  (kM  I  -S  )  -  0 

Jc  Jc  Jc 


i 
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Since  we  have  assumed  that  R,  t  0,  (4.9)  Implies  that 

Jc 

(4.10)  6  *  -  0 

and  thus  b  >  1.  Therefore  we  conclude  from  (4.5)  and  (4.10)  that 

(4.11)  (u.  *  AS  , 

~1  ^ 

c 

where  the  amplitude  A  is  to  be  determined. 

We  conclude  from  condition  (4.8)  with  j  =  2  and  (4.10)  and  (4.11)  that 
A  must  satisfy  the  equation 

(4.12)  or/i  A2  +  /  T)A  -  2r|'  -  0  , 

c  c  Jc 

where 

«  ■  ST  <FNH  +  2bn  !m  +  bn  +  2bn  bn  CNM>/(1  +  b„  \  >  • 

c  c  c  c  c  c  c 

it 

and  R  has  y.  *  \x  .  In  arriving  at  (4.12)  we  have  required  that  \l'  +  0 
■^c  c 

and  6"  +  0  which  by  (4.1)  imply  that  a  *  1  and  b  ■  2,  or 
(4.14)  M  “  Mc  +  ^lv  +  ,  5  =  v2 

There  are  two  real  roots  of  (4.12)  if 


+  ' 


(4.15)  Tf  >  - 


8  u  QfR. 
c  ], 
+2 
“l. 


Thus  there  are  two  inner  expansions 
r  N„ 


(4.16)  u> 


l«0J 


+  A* (71)5  (5)6%  +  0(5) 

Ml 

C 


I 


1 
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wlth 

(4.17)  M  -  +  T15*  +  0(5) 


If  orR+  >  0  Chen  these  expansions  are  defined  for  all  T|.  If  orR-!"  < 

3c  +<  +2  k 

there  are  two  expansions  for  T|  >  (-8fi  orR  /q.  )  and  two  expansions 

|C  c  c  l  ±| 

xl  9  L  I  4_9  X.  1  * 

7)  <  -  (-^lc0f Rj  /q‘  )*  At  710  =  (8 plaE* 
c  c  *  c 


/qp\ 

Jc 


dA 


d71 


0  then 
for 


21 


5.  The  Matching 

To  obtain  a  uniform  asymptotic  representation  of  the  solution  It  is 
necessary  to  connect  each  of  the  outer  expansions  to  the  appropriate 
inner  expansions  by  using  the  matching  conditions  of  the  method  of  matched 
asymptotic  expansions.  To  apply  these  conditions  we  first  express  the  outer 
expansion  in  terms  of  the  inner  variables  by  substituting  (4.14)  into  (3.1) 
and  expand  the  result  in  a  power  series  in  v: 

as  os 

(5.1)  u=  £  u.(g,#*(v))v2^  -  £  J&.(g)vk 

j-0  J  k-0 

The  matching  conditions  then  are 

(5.2)  lim  *  0  for  k  *  1,2,* •• 

hi-** 

For  the  outer  expansions  corresponding  to  the  uniform  sub-branches 
we  have 

+  ' 

-  \ 

(5.3) 

q,  c 

Jc 

and  corresponding  to  the  striation  sub-branches  we  have 


The  asymptotic  forms  of  the  inner  coefficients  iu^  are  found  from  the 
asymptotic  forms  of  the  amplitudes  A*  according  to  (4.11).  The  solution 
A*  of  (4.12)  is 
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1  -*  00 


/ 

R+' 

n 

-i- 

2 

U 

_ c 

+ 

% 

c 

lim  A  =  < 

*!)-*-• 

1 

+ 

^n 

_ c 

l  Mc 

a 

„-l 


of  >  0 


a  <  0 


+ 

In  obtaining  (5.6)  we  have  assumed  -  <  0  since  in  experiments  this  appears 

to  be  the  case. 

Hie  results  of  the  matching  are  summarized  in  Tables  I  and  II  where 

we  denote  the  uniform  and  striation  outer  expansions  for  pi  <  pi^  by  u^ 

•  4*  + 

and  respectively,  and  for  pi  >  by  u^  and  u^.  The  resulting 

composite  expansions  are  shown  in  Fig.  3. 

The  qualitative  features  observed  in  experiments  for  the  transition 

from  the  uniform  to  the  n  =»  1  striated  state  may  be  explained  by  Fig.  3a. 

The  fact  that  no  jump  in  the  state  of  the  plasma  is  observed  as  pi  is  varied 

rules  out  the  behavior  shown  by  Fig,  3b.  The  argon  contamination  and  non- 

uniform  heating  of  the  neon  plasma  are  therefore  sufficient  to  qualitatively 

explain  the  experimental  results. 

We  note  that  it  may  also  be  possible  for  a  weak,  steady  magnetic  field 
of  strength  B(z)  to  produce  smooth  transitions  between  states.  The  effect 
of  the  magnetic  field  is  to  reduce  the  radial  diffusion  of  electrons  in  the 
plasma  by  a  factor  which  depends  on  B(z).  The  term  in  F,  proportional  to  N, 
which  accounts  for  electron  radial  diffusion  losses  to  the  walls  of  the 
cylinder  will  therefore  have  a  coefficient  which  varies  with  z.  This 
effect  can  be  included  in  the  coefficient  in  (2.4a)  and  our  results 


remain  valid 
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Figure  Captions 


Fig.  1: 


Fig.  2: 


Fig.  3: 


The  contrast  between  theory  and  experiment  for 

stationary  striations.  The  amplitude  is  the 
n  ttz  + 

coefficient  of  sin  —  In  (1.2a)  with  /a  <  0. 

c 

The  curves  0^  and  0g  are  the  outer  expansions  corresponding 
to  the  uniform  and  striation  branches  of  the  bifurcation  problem, 
respectively.  The  dashed  curves  I  are  the  inner  expansions. 

The  amplitude  in  each  case  is  the  coefficient  of  sin  nc£. 

In  (a)  arR*  >  0  and  in  (b)  aR*  <  0. 

The  composite  solutions  of  the  perturbed  problem  for 

(a)  aR*  >  0  and  (b)  aR+  <  0.  The  noses  of  the  branches  in 

^c  *^c 

(b)  are  given  by  p  -  ±  T\Q6^  + 


•  •  • 


